On the origin of trigger-angle dependence of di-hadron 

correlations 



Wei-Liang Qian 1 , Rone Andrade 2 , Fernando Gardim 2 , Frederique Grassi 2 , Yogiro Hama 2 

1 Departamento de Fisica, Universidade Federal de Ouro Preto, MG, Brazil and 

2 Instituto de Fisica, Universidade de Sao Paulo, SP, Brazil 
(Dated: August 2, 2012) 

The STAR Collaboration reported measurements of di-hadron azimuthal correla- 
tion in medium-central Au+Au collisions at 200 A GeV, where the data are presented 
{3jq[ as a function of the trigger particle's azimuthal angle relative to the event plane <j) s . 

In particular, it is observed that the away-side correlation evolves from single- to 
double-peak structure with increasing <p s . In this work, we present the calculated 
correlations as functions of both (ft s and particle transverse momentum pj- , using 
the hydrodynamic code NeXSPheRIO. The results are found to be in reasonable 
agreement with the STAR data. We further argue that the above cp s dependence of 
the correlation structure can be understood in terms of one-tube model, as due to 
an interplay between the background elliptic flow caused by the initial state global 



geometry and the flow produced by fluctuations. 
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I. INTRODUCTION 



Di-hadron correlations in ultrarelativistic heavy-ion collisions provide valuable informa- 
tion on the properties of the created medium. The correlated hadron yields at intermediate 
and low pt, when expressed in terms of the pseudo-rapidity difference At] and azimuthal 
angular spacing A0, are strongly enhanced [l T3| compared to those at high pt 11. 12 . 
The structure in the near side of the trigger particle is usually referred to as "ridge" . It has 
a narrow A<fi located around zero and a long extension in At], and therefore it is tied to long 
range correlation in pseudo-rapidity. The away-side correlation broadens from peripheral to 
central collisions, and may exhibit double peak in A0 for certain centralities and particle pt 
ranges. The latter is usually called "shoulders". These structures, for the most part, can be 
successfully interpreted as consequence of collective flow jl3Ml 81] due to the hydrodynamical 
evolution of the system. 

Recently, efforts have been made to investigate the trigger-angle dependence of di-hadron 



correlations. STAR Collaboration reported measurements [19] of di-hadron azimuthal cor- 
relations as a function of the trigger particle's azimuthal angle relative to the event plane, 
<ps = \4>tr — ^ ep\ at different trigger and associated transverse momenta pr- The data are 
for 20-60% mid-central Au+Au collisions at 200 A GeV. In a more recent study [20] , the 
correlation was further separated into "jet" and "ridge" , where the ridge yields are obtained 
by considering hadron pairs with large \ Arj\. In this procedure, one assumes that the ridge is 
uniform in Arj while jet yields are not. Such assumption is quite reasonable when one takes 
into account the measured correlation at low px without a trigger particle j2l[. In their work, 
all correlated particles at \Arj\ > 0.7 are considered to be part of the ridge. It was observed 
that the correlations vary with <fr s in both the near and away side of the trigger particle. The 
"ridge" drops when the trigger particle goes from in-plane to out-of-plane and, moreover, 
the correlations in the away-side evolve from single- to double-peak with increasing <p s . 

Owing to the \Arj\ cut, it is quite probable that these data mainly reflect the properties of 
the medium. If this is the case, the main features of the observed trigger-angle dependence 
of di-hadron correlations should be reproduced by hydrodynamic simulations, as hydrody- 
namic models have been shown capable of reproducing many important characteristics in 



collective flow 
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The NeXSPheRIO code provides a good description of observed 

1 and two- 



collective flow[29l. l30l l36l . 137] . elliptic flow fluctuations 



pion interferometry [39] ; I n addition, it is known to reproduce the structures in di-hadron 
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long-range correlations [13|. In our previous studies on ridge 

mm 

1411 . we obtained some 

of the experimentally known properties such as: centrality dependence and pt dependence. 
It is therefore interesting to see whether the model is further able to reproduce the observed 
trigger-angle dependence of the data. Moreover, it is intriguing to identify the underlying 
physical origin behind the phenomenon and numerical simulations. We, therefore, propose 
an intuitive explanation on the mechanism of the trigger-angle dependence of the ridge 
structures based on hydrodynamics. This is the main purpose of the present study. 

The paper is organized as follows. In the next section, we carry out a hydrodynamic 
study on the trigger-angle dependence of di-hadron correlations by using NeXSPheRIO 
code. The calculations are done both with and without the pseudo-rapidity cut |A^| > 
0.7. Numerical results are presented for different angles of trigger particles and at different 
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associated-particle transverse momentum, and they are compared with STAR data [1 
We try to understand the origin of the observed features in section III, by making use of an 
analytic parametrization of one-tube model [13]. It is shown that the main features of the 
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experimentally observed S dependence can be obtained by the model, where one takes into 
account the interplay of the flow harmonics caused by a peripheral energetic tube and those 
from bulk background. 1 In our approach, the background modulation is evaluated both by 
using cumulant and ZYAM method, and very similar results are obtained. Section IV is 
devoted to discussions and conclusions. 



II. NUMERICAL RESULTS OF NEXSPHERIO 

Here, we present the numerical results on di-hadron azimuthal correlation, using the 
hydrodynamic code NeXSPheRIO. This code uses initial conditions (IC) provided by the 



event generator NeXuS [43|, |44j , solves the relativistic ideal hydrodynamic equations with 
SPheRIO code (45j. By generating many NeXuS events, and solving independently the 
equations of hydrodynamics for each of them, one takes into account the fluctuations of IC 
in event-by-event basis. At the end of the hydrodynamic evolution of each event, a Monte- 
Carlo generator is employed to achieve hadron emission, in Cooper-Frye prescription, and 
then hadron decay is considered. 

To evaluate di-hadron correlations, we generate 1200 NeXuS events in 20-60% centrality 
window for 200 AGeV Au-Au collisions. At the end of each event, Monte-Carlo generator 
is invoked for decoupling, from 300 times for more central to 500 times for more peripheral 
collisions. Here we emphasize that there is no free parameter in the present simulation, since 
the few existing ones have been fixed in earlier studies of r\ and px distributions [38]. To 
subtract the combinatorial background, we evaluate the two-particle cumulant. In order to 
make different events similar in characters, the whole centrality window is further divided 
equally into four smaller centrality classes, from 20-30% to 50-60%. Then one picks a trigger 
particle from one event and an associated particle from a different event to form a hadron 
pair. Averaging over all the pairs within the same sub-centrality class, one obtains the two 
particle cumulant. Background modulation is evaluated and subtraction is done within each 
sub-centrality class and then they are summed up together at the end. We calculated both 
cases: with and without the pseudo-rapidity cut. In the first case, all hadron pairs are 
included. Then in the latter case, only hadron pairs with |A?y| > 0.7 are considered as done 
in the STAR analysis. Note that in our approach, the IC constructed by "thermalizing" 
NeXus output do not explicitly involve jets, but they are not totally forgotten in the IC as 
they manifect themselves as high transverse fluid velocity in some localized region (see for 
instance Fig. 2 of Ref.j4o[). These regions with high transverse velocity, if not smeared out 
during hydrodynamic evolution, shall show up as a part of the near side "jet" peak in the 
resulting two particle correlaiton. Owing to its different physical origin, it is not correlated 
to the "ridge" and "shoulder" due to initial geometrical irregularities. The implementation 
of the pseudo-rapidity cut in our calculation should reduce the correlation due to such effect. 

The numerical results are shown in Fig.l and Fig.2 in solid lines. They are, respectively, 
compared with the STAR data pjl, 20] in filled circles and flow systematic uncertainties in 
histograms. From Fig.l and Fig.2, one sees that the data are reasonably reproduced by 
NeXSPheRIO code. Correlations decrease both in the near side and in the away side when 
cf) s increases. For high-pr triggers, this was thought to be related to the path length that 
the parton traverses 20]. Here, we see that such feature is also presented by intermediate- 



1 A preliminary report of this discussion has been presented by Y.H. at the ISMD2011 meeting. 42 1 
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and low-energy particles, and reproduced well in a hydro dynamic approach. The magnitude 
of correlations in Fig. 2 is globally smaller than those in Fig.l. This is because due to 
the Arj cut in the former case, the total yields of associated particles is reduced. If one 
only takes into consideration the overall shape of the correlations, one may notice that the 
hydrodynamic results in both plots are quite similar, Fig. 2 can be approximately obtained 
if one scales the plots in Fig.l by a factor of 0.6. The above results can be understood as a 
consequence of approximate Bjorken scaling in our hydrodynamic model. On one hand, the 
Arj cut effectively removes a portion of associated particles with selected pseudo-rapidity 
difference to the trigger particle. On the other hand, since the correlation is divided by the 
total trigger-particle number, the reduction of this does not affect the normalization. As 
expected, NeXSPheRIO results fit better for low momentum and the deviations increase at 
higher momentum. The calculated ridge structure in A0 varies, especially on the away side, 
with trigger direction. For in-plane trigger, simulation results exhibit a one-peak structure, 
which is broader in comparison with the near-side peak. The away-side structure changes 
continuously from one peak at in-plane direction (<f) s = 0) to double peak at out-of-plane 
direction {(f) s = tt/2). This characteristic is manifested particularly for larger transverse 
momentum of the associated particle. All these features are in good agreement with the 
STAR measurements. 



III. THE ORIGIN OF TRIGGER- ANGLE DEPENDENCE OF RIDGE 

STRUCTURE 

As pointed out, the motivation to introduce the \Ar]\ cut was to separate "ridge" from 
"jet", therefore the measured "ridge" correlations are expected to reflect the properties of 
the medium. In the previous section, we showed that hydrodynamic simulations are able 
to reproduce the main features observed experimentally. Now, how this effect is produced? 
In a previous study [46j], we tried to clarify the origin of the effect by using the one-tube 
model fl4| adapted to non-central collisions. Full details of the model can be found in 
[Hi fl8i l40l l4li 47]. Here we only outline the main thoughts. In one-tube model, or more 



generally peripheral-tube model (allowing more than one tube), one treats the initial energy 
profile on a transverse plane as superposition of peripheral high energy tubes on top of 
the background, which can be thought as the average distribution. The problem is further 
simplified by studying the transverse hydrodynamic expansion of a system consisting of 
only one peripheral tube with the assumption of longitudinal invariance. As the high energy 
tube deflects the collective flow generated by the background, the resulting single particle 
azimuthal distribution naturally possesses two peaks. They eventually give rise to the desired 
"ridge" and "shoulder" structures. In ref. [46], we took as background the average energy- 
density distribution obtained with NeXSPheRIO, which has an elliptical shape for non- 
central collisions. A peripheral tube sits on top of the background and its azimuthal position 
varies from event to event. As shown there, the di-hadron correlation obtained from such a 
simple IC configuration does reproduce the main features of the data. 

Though the IC of the problem was greatly simplified through such an approach, the 
underlying physical mechanism of the obtained features were still not clear. In order to 
identify it more transparently, we further simplify it using an approximate analytical model. 
The derivation of the results relies on the following three hypotheses: 



• The collective flow consists of contributions from the background and those induced 



5 



by a peripheral tube. 

• A small portion of the flow is generated due to the interaction between background 
and the peripheral tube and, therefore, the flow produced in this process is correlated 
to the tube. 

• Event by event multiplicity fluctuations are further considered as a correction that sits 
atop of the above collective flow of the system. 

Let us comment briefly on these hypotheses. Based on the idea of one-tube model, a small 
portion of the background flow is deflected by the peripheral tube; extra Fourier components 
of the flow are generated by this process. The event plane of these extra flow harmonics 
are consequently correlated to the location of the tube, as stated in the second hypothesis. 
Since the contribution from the tube is small, we will treat it perturbatively, considering 
the resultant flow a superposition of the background flow and the one produced in the tube- 
background interaction as described above. We believe that, at least qualitative behavior 
of the results will remain valid also in more realistic case. Here, we are considering just one 
tube as fluctuation. However, in the limit of small perturbations, it is quite straightforward 
to generalize our results to the case of N tubes. There are many possible sources for 
fluctuations, such as: flow fluctuations, multiplicity fluctuations etc, we will only consider 
multiplicity fluctuations in this simple model as assumed in the third hypothesis. As it will 
be shown below, this turn out to be enough to derive the observed feature in di-hadron 
correlations. 

Using the hypotheses stated above, we write down the one-particle distribution as a sum 
of two terms: the distribution of the background and that of the tube. 

dN dN bgd dN tube 



(0) + ^ (0)0t); (1) 



where 



dNbgd ( , s _ N b 



(0) = -1(1 +^008(20)), (2) 

tup Z7T 

^T-Mt) = 7T E 2^cos(r#-0 t ]) (3) 

Since the background is dominated by the elliptic flow in non-central collisions, as observed 
experimentally, in Eq.(J2]) we consider the most simple case, by parametrizing it with the el- 
liptic flow parameter v\ and the overall multiplicity, denoted by N b . As for the contributions 
from the tube, for simplicity, we assume they are independent of its angular position <p t and 
take into account the minimal number of Fourier components to reproduce the two-particle 
correlation due to the sole existence of a peripheral tube in an isotropic background (see 



the plots of Fig.2 of ref. 18|, for instance). Therefore, only two essential components v\ and 
v\ are retained in Eq.Q. We note here that the overall triangular flow in our approach is 
generated only by the tube, i.e., v 3 = v\ and so its symmetry axis is correlated to the tube 
location <j) t . The azimuthal angle <p of the emitted hadron and the position of the tube <j>t 
are measured with respect to the event plane \&2 of the system. Since the flow components 
from the background are much bigger than those generated by the tube, as discussed below, 
\&2 is essentially determined by the elliptic flow of the background v\. 
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The di-hadron correlation is given by 



proper I J AT \ m "i- xe d 



>Hw w ) -(w ( «; • <4) 

where <p s is the trigger angle (<f) s = for in-plane and (f) s = it/ 2 for out-of-plane trigger). In 
one-tube model, 



J at \ proper 
u ' 1 Vpair 



f dd> t „. , .dN , , , .dN . a i i \ 

where f(4>t) is the distribution function of the tube. We will take f(4>t) = 1, for simplicity. 

The combinatorial background (dN pa i r / dA(j)) mixed can be calculated by using either cu- 
mulant or ZYAM method |48j. As shown below, both methods lend very similar conclusions 
in our model. Here, we first carry out the calculation using cumulant, which gives 

/dN pair \ ™« f d<j> t fU v f W t f( ,,,dN dN 

Notice that, in the averaging procedure above, integrations both over <ft t and <fi' t are required 
in the mixed events, whereas only one integration over <p t is enough for proper events. This 
will make an important difference between two terms in the subtraction of Eq.(j3]). 

Using our simplified parametrization, Eqs. (tUlHl) and, by averaging over events, one obtains 

dN ■ \ proper < N 2 > 

x b a + 2^cos(20 s ))(l + 2^cos(2(A0 + s ))) 



dA(f) I (2tt) 2 

■Nt 
2tt 



+ (§) 2 E 2 «) 2cos (^) ( ? ) 

n=2,3 



and 



J at \ mixed(cmlt) at 2 

= -^^(l + 2^cos(20 s )(l + 2^cos(2(A0 + ^))). (8) 

Because of random distribution, contributions from peripheral tube are cancelled out upon 
averaging in the mixed-event correlation. Observe the difference between the factors mul- 
tiplying the background terms of the proper- and mixed-event correlation. By subtracting 
Eq.(JS]) from Eq.([7j), the resultant correlation is 

J AT \ {emit) N 2 ^ _ ^ AT ^2 

A b (l + 2^cos(20 s ))(l + 2^cos(2(A0 + 0,))) 



dA(f) v ^ s 7 (2tt) 

.Nt 
'2tt 



+ (^) 2 E 2 «) 2cos (^)- ( 9 ) 

n=2,3 



So, one sees that, as the multiplicity fluctuates, the background elliptic flow does contribute 
to the correlation. Now, by taking <p s = in Eq.Q, the correlation for the in-plane trigger 
is finally given as 

dN pair \ (c "" ,> <A?>-<iV t >* , , . 



/ )2 ;i + 2„;)(l + 2^cos(2A*)) 

' / in— plane V / 

+ (§) 2 E 2 (^) 2cos (^)- ( 10 ) 

n=2,3 
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We note that, due to the summation of the two terms concerning v\ and t^, the away-side 
peak becomes broader than the near-side one, as shown below in Figj3l 
Similarly, the out-of-plane correlation is obtained by putting (fi s = n/2 as 

dA-V-" _ < N £ > - < _ _ Hms(2m) 



dA<j) I out _ of _ vlane (2tt) 

.Nt 
'2tt 



+ (^) 2 E 2 «) 2cos (^)- ( n ) 



n=2,3 

One sees that, because of the shift in the trigger angle (fi s (0 — > 7r/2), the cosine depen- 
dence of the background contribution gives an opposite sign, as compared to the in-plane 
correlation. This negative sign leads to the following consequences. Firstly, there is a re- 
duction in the amplitude of out-of-plane correlation both on the near side and on the away 
side. More importantly, it naturally gives rise to the observed double peak structure on the 
away side. Therefore, despite its simplicity, the above analytic model reproduces the main 
characteristics of the observed data. The overall correlation is found to decrease meanwhile 
away-side correlation evolves from a broad single peak to a double peak as 4> s increases. The 
correlations in Fig. 3 is plotted with the following parameters 

< N? >= 0.45 
4 = 4 = 0.1 
< N% > - < N b > 2 = 0.022 

v\ = 0.25 (12) 



We note that both the correlated yields and the flow harmonics are actually dependent on 
the specific choice of the pt interval of the observed hadrons as shown in Fig.l and Fig.2. 
Since the transverse momentum dependence has not been explicitly taken into account in 
this simple model, the multiplicities in the above parameter set are only determined up to an 
overall normalization factor, and the flow coefficients are chosen to reproduce the qualitative 
behavior of the trigger-angle dependence of di-hadron correlation as shown by data. 

Now we will show that very similar results will be again obtained, if one evaluates the 
combinatorial mixed event contribution using ZYAM method. The spirit of ZYAM method 
is to first estimate the form of resultant correlation solely due to the average background 
collective flow and, then, the evaluated correlation is rescaled by a factor B, the latter is 
determined by assuming zero signal at minimum of the subtracted correlation. Di-hadron 
correlation for the background flow is given by 

dN pair \™ xed(ZYAM) [d<i> x ,, , sdN bgd dN bgd 



dA0 



b{*.) J ^-^ dJ W^ dJ W^ + ^- (13) 



In the STAR analyses, both t>2 and v\ have been taken into account for the background. In 
our simplified model, however, the average background flow contains, besides the radial one, 
only the elliptic flow v\. Therefore, a straightforward calculation gives 

/ J AT \mixed(ZYAM) jy 2 

{-^) = B &')^£~0- + 2^cos(20 s )(l + 2v\ cos(A0 + <f> s )). (14) 
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We remark that, since in ZYAM method fluctuations are not explicitly considered, the 
multiplicity of background distribution, as given by Eq. ([2|) , is evidently average multiplicity 
< Nb >. By combining this term with the proper correlation, given by Eq.((7]), the resultant 
correlation reads 

dN pair ^.\ {ZYAM) < N 2 b > -B{<j> 8 ) < N b 



dA(f) vrs 7 (2?r) 

'2tt 



[1 + 1v\ cos(20 a ))(l + 2v\ cos(2(A0 + S ))) 

+ (^) 2 E 2 ^ 2cos (^) ( 15 ) 



n=2,3 

The consistency of this expression with what is used in the STAR analyses will be discussed 
below, in the next section. Note that the normalization factor B(<f) s ) is a function of trigger 
angle. It is fixed to give zero yield at the minimum of the subtracted correlation, namely 
(dN pair /dA^ s )) {ZYAM) = at the minimum. Because the correlation in Eq.f ll5p is posi- 
tively defined, and the second term can be positive or negative, the coefficient of the first 
term < N% > —B((p s ) < N b > 2 must be positive. One sees clearly that the above expression 
is almost identical to the cumulant result, Eq.Q, so will cause similar trigger-angle depen- 
dence. This can also be seen from the plots of correlations obtained by adopting the same 
parameters as in Eq.f fT2~]) and additionally 

< N% >= 100, (16) 

where the value of < N% > is choosen to be larger compared to its fluctuation. The scale 
factor B(<p s ) is subsequently fixed by the minimum condition as 

B((j) s = 0) = 1.000053 

B(<p s = tt/2) = 1 (17) 

The in-plane correlation plot is shown in Fig.4, which is very close to the one in Fig.3. The 
one-tube contribution and out-of-plane correlation are not plotted, since they are exactly 
the same as those shown in Fig.3. 



IV. DISCUSSIONS AND CONCLUSIONS 



Here, we first show that the expression of di-hadron correlation of background flow in 
Eq. (1141) is in agreement with that obtained in Ref . [49( , which is employed in STAR analyses 
20] . The only difference is that we have only taken into account the second-order har- 
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monic, and the reason for not including any higher-order flow components in our calculation 
is simply because we wanted to transparently show the mechanism of in-plane/out-of-plane 



effect by using a model as simple as possible. In Ref. 49], it was shown that 



dNpair 

dA(f) 



B (R) 



l + 2j2 V n )v n R) COs(nA0) 



n=l 



(18) 



where is the associated particle's n-th harmonic, Vn is the average n-th harmonic 
of the trigger particles, and B^ R \ the background normalization, denotes the integrated 
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inclusive pair yield 

sin (fee) 



B (R) = 51+ Yl 2v k ] cos ( k <f>s 



fe=2,4,6,. 



kc 



,, D , V^n + ^n.even^n + J2k=2A,6,... [ V k+n + %-ral J 

4 ' } = \r, , 19 

1 4- V 2v ( ) Tu 

m ., , . sin(A;c) . . T ,. 

T fc = cos(A;0 s )— ^ (cos(A;A$)> , 

where 2c is the angular width where a trigger is located. In our approach, the size of the slice 
is taken to be infinitely small (c — > 0), and perfect event plane resolution ((cos(A;A\l/)) = 1) is 
assumed. Take, for instance, the in-plane correlation by substituting S = 0, and considering 
terms up to the second order, one obtains 



BW = B (l + Ivf 



vf R) = 1 

which is readily shown to be consistent with Eq.( TH| . In fact, it is intuitive to understand 
since, the one-particle distribution of the trigger in this case is a 5 function which peaks at 
S = 0. 

In our approach, the trigger-angle dependence of di-hadron correlation is understood as 
due to the interplay between the elliptic flow caused by the initial almond deformation of the 
whole system and flow produced by fluctuations. The contributions due to fluctuations are 
expressed in terms of a high-energy-density tube and the flow deflected by it. However, the 



generic correlation due to the tube is preserved even after the background subtraction [18 
by this simple model we show explicitly that the result does not depend on either cumulant 
or ZYAM method. This is because the form of combinatorial background is determined by 
the average flow harmonics, as shown in Eqs.flU]) and ffl5l) . Though in this approach, only 
elliptic flow is considered for simplicity, it is straightforward to extend the result here to 
a more general case. Due to multiplicity fluctuations or due to the procedure in ZYAM, 
the background flow may also contribute to the subtracted correlation. Since the back- 
ground modulation is shifted, changing the phase, when the trigger particle moves from 
in-plane to out-of-plane as seen in EqsflS]) and (TH|) . the summation of the contributions of 
the background and that of the tube give rise to the desired trigger-angle dependence. 

In the one-tube model, a part of the flow is caused by the peripheral energetic tube. Since 
the tube deflects the global flow of the background, the event planes of such flow components 
are correlated to the localization of the tube as expressed in Eq.(j3]). Particularly, it also 
contains the second harmonic v\. Though it is present in the proper two particle correlation 
(Eq.(E])), it is not considered in Eq. (fl4|) when evaluating combinatorial background correla- 
tion. The reason is two- fold. Firstly, to calculate the average t>2 of a given event, one must 
use multiplicity as weight, in our model, the multiplicity of background N b is assumed to be 
much bigger than that of the tube N t . (The parameter < N£ >= 100 can be freely changed 
to a much bigger number.) Moreover, since t varies from event to event, the contribution 
to t>2 from v\ is positive at <p t = and negative at t = tt/2. When averaging over different 
events, most contributions cancel each other at different <p t values. As it happens, v\ con- 
tributes to the subtracted correlation while it does not manifest itself in average background 
flow. This is an important feature of the present approach. 
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It is interesting to note that the two particle correlation has also been studied using 



Fourier expansion in[17|, |50 



dN v proper ^ 

^f^s)) = ^(1 + 2V 2A cos(2A0) + 2V 3 acos(3A0)) + --- (20) 
For comparison, we rewrite Eq. ((Tj) in terms of V n A as follows 

V ^ = Tm\ hJSb Hf + cos(2^)«S (21) 

(N 2 ) (l + 2v°cos(20J) 

N 2 

V3A = (JV*> (l + 2v>cos(20 fl )) • (22) 

One sees that the background elliptic flow v\ dominates V2A for both in-plane and out-of- 
plane directions, while V3A is determined by the triangular flow v\ produced by the tube. 
Due to the factor cos(20 s ), the second term of Eq. (j2T]) changes sign when the trigger angle 
goes from (p s = to 4> s = it/2. Dominated by this term, V2A decreases with <f) s , and it 
intersects V2A — at around <p s = tt/4. Since the first term in Eq.( l21j) is positive definite, 
the integral of V2A with respect to <p s is positive. These features are in good agreement 
with the data analysis (see Fig.l of ref.(l7j|). On the other hand, the axis of triangularity is 
determined by the position of the tube ^3 = <p t . Since we have assumed uniform distribution 
/(0t) = 1 in our calculation, the event plane of triangularity is totally uncorrelated with the 



event plane ^2, as generally understood [151 . 15 11 . and consequently, the contribution from 



triangular flow should not depend much on the event-plane angle. This is indeed shown in 
the above expression Eq (f22]) . V 3 a barely depends on <p s , if anything, it slightly increases 
with increasing <p s . This characteristic is also found in the data 171 . 



In conclusion, the NeXSPheRIO code gives correct qualitative behavior of the in- 
plane/out-of-plane effect. Physically, we understand that this effect appears because, besides 
the contribution coming from the peripheral tube, additional contribution from the back- 
ground flow has to be considered. The latter is back-to-back (peaks at A0 = 0, 7r) in the 
case of in-plane triggers (<f> 8 = 0) and rotated by 7r/2 (peaks at A0 = — 7r/2,7r/2) in the 
case of out-of-plane triggers (<f> s = vr/2). A simplified analytical model is proposed, and it 
is shown to successfully reproduce the observed features. 
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FIG. 1: The subtracted di-hadron correlations as a function of Acp for different <p s = (fitrig — 4>EP 
and p^ soc with 3 < p^ 19 < AGeV in 20 - 60% Au+Au collisions. The 4> s range increases from 0-15° 
(left column) to 75-90° (right column); the p9f soc range increases from 0.15-0.5 GeV (top row) to 
1.5-2 GeV (bottom row). NeXSPheRIO results in solid curves, are compared with STAR data in 
filled circles [ijj. The histograms indicate the systematic uncertainties from flow subtraction. 
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FIG. 2: The subtracted di-hadron correlations as a function of Acp for different <p s = (fitrig ~ 4>EP 
and pf soc with 3 < p% i(J < AGeV and |A??| > 0.7 in 20 - 60% Au+Au collisions. The <j) s range 
increases from 0-15° (left column) to 75-90° (right column); thep^ soc range increases from 0.15-0.5 
GeV (top row) to 2-3 GeV (bottom row). NeXSPheRIO results in solid curves are compared with 
STAR data in filled circles [2oT |. The grey histograms indicate the systematic uncertainties from 
flow subtraction. 
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FIG. 3: Plots of di-hadron correlations calculated by cumulant method. From the left to the right: 
(i) the peripheral-tube contribution; (ii) the one from the background (dashed line) and the resul- 
tant correlation (solid line) for in-plane triggers, as given by Eq. (|10p ; and (hi) the corresponding 
ones for the out-of-plane triggers, Eq. (fTT|) , 
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FIG. 4: Plot of in-plane di-hadron correlation by using ZYAM as given by Eq. (]15p , The one-tube 
contribution and out-of-plane correlation are exactly the same as the ones shown in Fig. 3. 



